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The production of K*(892)° and 0(1020) in pp collisions at ^/s=l TeV was measured by the ALICE 
experiment at the LHC. The yields and the transverse momentum spectra d 2 N/dydpj at midrapidity 
|y|< 0.5 in the range 0<p T <6 GeV/c for K*(892)° and 0.4<p T <6 GeV/c for 0(1020) are reported 
and compared to model predictions. Using the yield of pions, kaons, and Q. baryons measured previ- 
ously by ALICE at ^=7 TeV, the ratios K7K", 0/K*, 0/K", <j)/n~, and (il + £2)/0 are presented. 
The values of the K*/K~, 0/K* and 0/K~ ratios are similar to those found at lower centre-of-mass 
energies. In contrast, the <p/K~ ratio, which has been observed to increase with energy, seems to 
saturate above 200 GeV The (£2 + i2)/0 ratio in the pj range 1-5 GeV/c is found to be in good 
agreement with the prediction of the HIJING/BB v2.0 model with a strong colour field. 
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1 Introduction 

The study of resonance production plays an important role in both elementary and heavy ion collisions. 
In pp and e + e~ collisions it contributes to the understanding of hadron production [T] [2] as the decay 
products of resonances represent a large fraction of the final state particles. In addition, it provides 
a reference for tuning event generators inspired by Quantum Chromodynamics (QCD). In heavy ion 
collisions, resonances are a sensitive probe of the dynamical evolution of the fireball. Due to their short 
lifetime (a few fm/c) a significant fraction of resonances decay inside the hot and dense medium and 
their hadronic daughters interact with the medium during the fireball expansion (3] SI 13 . 

The 0(1020), which is the lightest vector meson composed only of sea quarks, provides a probe for the 
study of the strangeness production. In pp collisions, ss pair production was found to be significantly 
suppressed in comparison to uu and cW-pair BIT]]. Another useful probe of strangeness production is 
the K*(892)°, which is a vector meson with a mass similar to the 0, but differing by one unit of the 
strangeness quantum number. The (Q. + Q.)/(j) ratio has been suggested [ 8 ] as a probe of the colour field 
strength, which in microscopic models influences the relative yield of strange with respect to non-strange 
particles. 

We present the first measurement of the differential (d 2 N/dydpj) and /^-integrated (dN/dy) yields of 
the K* Q and 0(1020) mesons at midrapidity (|y|<0.5) in pp collisions at ^=1 TeV. The data analysis 
was carried out for K* (0) on a sample of 80 (60) million minimum bias pp collisions collected by the 
ALICE experiment. The resonances were identified via their main decay channel K* — > n ± +K T and 
— > K++KT. Tracks were reconstructed by the main ALICE tracking devices, the Time Projection 
Chamber (TPC) and the Inner Tracking System (ITS). The TPC and Time of Flight (TOF) detectors 
were used to identify pions and kaons. The measured spectra are compared to two QCD-based event 
generators, PHOJET El and PYTHIA iflOl . 

The ratios K*/K~, 0/K*, 0/K~, and 0/7T - are computed using the yield of pions and kaons measured ifTTTl 
with the ALICE detector in pp collisions at 7 TeV. These ratios are compared with measurements at 
lower collision energies. The (Q. + Q.) /0 ratio has been calculated as a function of transverse momentum 
using the Q. and Q. yield measured at 7 TeV [12]; this ratio is then compared to the predictions of the 
HIJING/BB v2.0 model with a Strong Colour Field (SCF) (H and to PYTHIA-Perugia 2011 Ifl4ll . 

The article is organized as follows: Section 2 gives details about the detectors relevant for this analysis, 
Section 3 describes the criteria used for event selection, Section 4 gives an overview of the analysis, 
Section 5 presents the results and Section 6 the conclusions. 

2 Experimental set-up 

A full description of the ALICE detector can be found in |[I31[T6i l. For the analyses described in this paper, 
the ITS, the TPC, and the TOF detectors were used. These detectors are set inside a large solenoidal 
magnet providing a magnetic field B=0.5 T, and have a common pseudorapidity coverage of |r/| < 0.9. 
Two forward scintillator hodoscopes (VZERO) placed along the beam direction at -0.9 m and 3.3 m on 
either side of the interaction point, cover the pseudorapidity regions —3.7 < tj < — 1.7 and 2.8 < tj < 5.1. 
These are used for triggering and for rejecting beam-gas interactions. 

2.1 The Inner Tracking System 

The ITS Ifl6l is the innermost ALICE detector, located between 3.9 and 43 cm radial distance from the 
beam axis. It is made of six cylindrical layers of silicon detectors (two layers of pixels, two of silicon 
drift, and two of silicon strips), with a total material budget of 7.66 % of the radiation length Xq. It 

1 We denote by K* the average of K*(892)° and K* (892)°. 
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provides high-resolution space points close to the interaction vertex, thus improving momentum and 
angular resolution of the tracks reconstructed in the TPC. 

The two innermost ITS layers constitute the Silicon Pixel Detector (SPD), which has a high granularity 
of about 9.8 million pixel cells, each with a size of 50x425 jUm 2 . These layers are located at radii of 3.9 
and 7.6 cm with pseudorapidity coverages of |t] | < 2.0 and \t] \ < 1 .4, respectively. The detector provides 
a position resolution of 12 in the r(p direction and about 100 in the direction along the beam 
axis. 

2.2 The Time Projection Chamber 

The TPC [ 17] is the main ALICE tracking device. It is a large- volume, high-granularity, cylindrical drift 
detector which has a length of 5.1 m and inner and outer radii of 0.85 and 2.47 m, respectively. It covers 
the pseudorapidity range | r\ \ < 0.9 with a full azimuthal acceptance. The drift volume is filled with 90 m 3 
of Ne/C02/N2. The maximum drift time is 94 /is. A total of 72 multi-wire proportional chambers with 
cathode pad readout instrument the two end plates, which are segmented into 18 sectors and include a 
total of over 550,000 readout pads. The ionization electrons drift for up to 2.5 m and are measured on 159 
pad rows. The momentum resolution of the TPC is in the range 1-7% for pions with 1</;>t<10 GeV/c. 
The ALICE TPC ReadOut (ALTRO) chip, employing a 10 bit ADC at 10 MHz sampling rate and digital 
filtering circuits, allows for precise position and linear energy loss measurements with a gas gain of the 
order of 10 4 . The material budget of the TPC near T] = amounts to about 4.1% of Xq. 

The position resolution in the r(p direction varies between 1100 /im and 800 /im going from the inner to 
the outer radius, whereas the resolution along the beam axis varies between 1250 /im and 1100 /im. 

2.3 The Time Of Flight detector 

The ALICE TOF ifTSl [T9l is a cylindrical assembly of Multi-gap Resistive Plate Chambers (MRPC) with 
an inner radius of 370 cm and an outer radius of 399 cm. It has a pseudorapidity coverage of \r] \ < 0.9 
and full azimuthal acceptance, except for the region 260° < (p < 320° at |t] | < 0.14 where a gap was left 
in order to reduce the amount of material in front of the Photon Spectrometer (PHOS). The elementary 
unit of the TOF system is a 10-gap double-stack MRPC strip 122 cm long and 13 cm wide, with an 
active area of 120x7.4 cm 2 subdivided into two rows of 48 pads of 3.5x2.5 cm 2 each. The length of 
the TOF barrel active region is 741 cm. It has about 153,000 readout channels and an average thickness 
of 25-30% of Xo, depending on the detector zone. For pp collisions, such a segmentation leads to an 
occupancy below 0.02%. The front-end electronics are designed to comply with the basic characteristics 
of the MRPC detector, i.e. very fast differential signals from the anode and the cathode readout: the 
resulting intrinsic time resolution of the detector and electronics was measured to be smaller than 50 ps. 

3 Data collection and event selection 

Data used for this analysis were collected in 2010 using a magnetic field of B=0.5 T with both field 
polarities. The minimum bias trigger required a single hit in the SPD detector or in one of the two 
VZERO counters, i.e. at least one charged particle anywhere in the ~8 units of pseudorapidity covered 
by these detectors. In addition, a coincidence was required with signals from two beam pick-up coun- 
ters, one on each side of the interaction region, indicating the passage of proton bunches. The trigger 
selection efficiency for inelastic collisions was estimated to be 85.2% with a +6.2% and -3% relative 
uncertainty 1201 . During the data-taking period, the luminosity at the ALICE interaction point was kept 
in the range 0.6 — 1.2 x 10 29 cm _2 s _1 . Runs with a mean pile-up probability per event larger than 5% 
were excluded from the analysis. 

Beam-induced background was reduced to a negligible level (< 0.01%) with the help of the timing 
information of the VZERO counters and by a cut on the position of the primary vertex reconstructed 
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by the SPD Gil . Accepted events were required to have a reconstructed primary vertex. Its position 
can be computed either using the tracks reconstructed by TPC and ITS, or using the "tracklets" obtained 
connecting reconstructed clusters in both SPD layers. If possible, the first method is used. First, for 
each event a three dimensional reconstruction of the primary vertex was attempted with either a Kalman 
filter, using reconstructed tracks as input, or by a minimization of the squared distances between all 
the extrapolated tracklets. Otherwise only the z position of the primary vertex was reconstructed by 
correlating the z coordinates of the SPD space points, while for x and y the average position of the beam 
in the transverse plane was taken. The primary vertex reconstrunction efficiency, calculated via Monte 
Carlo simulation, approaches unity in events with a K* or a <p produced in the central rapidity region. In 
order to minimize acceptance and efficiency biases for tracks at the edge of the TPC detection volume, 
events were accepted only when their primary vertex was within ±10 cm from the geometrical centre of 
the ALICE barrel. 

4 Data analysis 

4.1 Track selection 

Global tracking in ALICE is performed using ITS and TPC clusters. It is based on a Kalman filter al- 
gorithm which takes into account both multiple scattering and energy loss along the path as described in 
detail in 11221 . The Distance of Closest Approach (DC A) to the primary vertex is used to discriminate 
between primary and secondary particles. Primary charged particles are those produced directly in the 
interaction and all decay products from particles with a proper decay length cz <1 cm; secondary parti- 
cles include those from the weak decay of strange hadrons and from interactions in the detector material. 
Several cuts were applied to achieve a high track quality in the analyzed sample. Tracks were required 
to have at least 70 reconstructed clusters in the TPC out of the maximum 159 available. This ensured 
a high efficiency and good d£Vdx resolution, keeping the contamination from secondary and fake tracks 
small. 

In order to improve the global resolution, tracks were accepted only in the range \ < 0.8 (i.e. well 
within the TPC acceptance) and with pj> 0.15 GeV/c. In order to reduce secondary particles, tracks 
were required to have at least one hit in one of the two innermost backing detectors (SPD) and to have 
a DCA to the primary vertex less than 2 cm along the beam direction. The DCA in the transverse plane 
was required to be smaller than 7 <7dca(/?t), where cTdcaQ?t) = (0.0026 + 0.0050 GeV/c -pj _1 ) cm. 

4.2 Particle identification 

Identification of pions and kaons is performed using the measurements of the TPC and the TOE For the 
TPC, the particle is identified based on the energy it deposits in the drift gas, compared with the expected 
value computed using a parameterized Bethe-Bloch function Il23l 1241 . Figure Q] shows the TPC signal 
versus track momentum computed at the point the particle enters the detector, and the curves represent 
the Bethe-Bloch functions for each mass hypothesis. The TPC calibration parameters have mostly been 
determined and tested via the analysis of cosmic rays; the chamber gain has been measured using the 
decay of radioactive 83 Kr gas released into the TPC volume ifTTl . 

A truncated-mean procedure is used to determine d£Vdx, with only 60% of the points kept. The d£Vdx res- 
olution Otpc is about 5% for tracks with 159 clusters and about 6.5% when averaged over all recon- 
structed tracks. The relevant value of Otpc is estimated for each track taking into account the actual 
number of clusters used IfTTl . 

The TPC d£Vdx measurement allows pions to be separated from kaons for momenta up to p ~ 0.7 GeV/c, 
while the proton/antiproton band starts to overlap with the pion/kaon band at p «1 GeV/c. As can be 
observed in Fig. [TJ the electron/positron dE/dx band crosses the other bands at various momenta. This 
contamination in identified pions and kaons can be drastically reduced using information from the TOE 
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Fig. 1: (Colour online) Specific ionization energy loss dE/dx vs. momentum for tracks measured with the ALICE 
TPC The solid lines are parametrizations of the Bethe-Bloch function [23 1. 

Particles are identified in the TPC via the difference between the measured energy loss and the value 
expected for different mass hypotheses. The cut on this difference, normalized to the resolution Gjpc, 
is optimized for each analysis and depends in general on the signal-to-background ratio and on the 
transverse momentum. 

Figure |2] shows the correlation between particle momentum and their velocity /3 = L/ct, where L is the 
total integrated path length and t is the time of flight measured by the TOF detector. For the analyses 
described in this paper the start time of the collision is estimated using the particle arrival times at the 
TOF or the averaged collision time observed in the fill. The bands corresponding to pions, kaons, protons 
and deuterons are clearly visible. 

Particles are identified in the TOF by comparing the measured time of flight to the expected time for 
a given particle species. The cut is expressed in units of the estimated resolution Gtof for each track, 
which has a mean value of 160 ps. The TOF allows pions and kaons to be unambigously identifed up to 
p ~ 1.5 GeV/c. The two mesons can be distinguished from (anti)protons up to p ~ 2.5 GeV/c. 

Considering the high multiplicities reached in pp collisions at LHC energies, good particle identification 
is important to reduce combinatorial background as well as correlated background from misidentified 
resonance decays. The <p analysis requires only primary kaons to be selected and cuts were kept loose 
in order to maximize the efficiency. The cut for particle identification in the TPC was set to 3gtpc 
(5c7tpc) for tracks with p larger (smaller) than 0.35 GeV/c. When a TOF signal is present, a particle 
identification cut of 30tof is also applied. For the K* analysis, both pions and kaons are identified. Two 
different strategies were followed. For tracks with TOF signals, a TPC d£Vdx cut of 5c7tpc was applied 
and a TOF cut of 3ctof (2otof) was a PP ue d for tracks with momenta below (above) 1.5 GeV/c. For 
tracks without a TOF signal the kaon momentum was required to be below 0.7 GeV/cand 5c7tpc, 3gtpc, 
and 2gtpc cuts were used for p <0.35 GeV/c, 0.35< p < 0.5 GeV/c, and p > 0.5 GeV/c, respectively. 
This more restrictive cut on kaons was used to reduce the correlated background originating from p 
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Fig. 2: (Colour online) Velocity j3 of particles measured by TOF vs. momentum, 
decays in which a pion is misindentified. 
4.3 Raw yield extraction and background estimation 

The uncorrected background was estimated using two different techniques: like-sign and event mixing. 
In the like-sign method invariant mass distributions of like-sign Ktt or KK combinations (for K* and 
0, respectively) from the same event were constructed. In the event mixing method the shape of the 
uncorrected background was estimated from the invariant mass distribution of unlike-sign Kn or KK 
combinations from different events. To avoid mismatch due to different acceptances and to assure a 
similar event structure, tracks from events with similar vertex positions z (Az < 1 cm) and track mul- 
tiplicities n (An < 10) were mixed. To reduce statistical uncertainties each event was mixed with 10 
other events. The mixed-event distribution was then normalized in the mass region 1.08 < M <1.2 (1.04 
< M <1.07) GeV/c 2 for K* (0), and subtracted in each pj bin. The uncertainty in the normalization was 
estimated by varying the normalization region and is included in the quoted systematic uncertainty for 
signal extraction. After background subtraction a residual background remains. This is due in part to 
an imperfect description of the combinatorial background but mainly caused by a real correlated back- 
ground. The latter can arise from correlated 7rK or KK pairs or from misidentified particle decays (for 
example K*° for (p, or and p for K*, or from underlying jet event structure). 

The total /^-integrated number of reconstructed mesons after background subtraction was about 1.8x 10 6 
for the K* and 2.3 xlO 5 for the 0. For the K* the signal-to-background ratio varied from 0.08 at 
/?t=0.05 GeV/c to 0.2 at pj= 5.5 GeV/c. The significance (S/y/S+B) was about 34 in the pj bins 
at both 0.05 and 5.5 GeV/c and reached a maximum of about 127 at 1 GeV/c. For the (j) the signal-to- 
background ratio varied from 2.8 to 1.6 between /?t=0.45 and /?t=5.5 GeV/c, with a minimum of 0.5 at 
1.6 GeV/c; the significance was about 30 in the pj bins at both 0.45 and 5.5 GeV/c with a maximum of 
90 at 1 GeV/c. 

The raw yield of K*(892)° and its antiparticle was extracted in 22 pj bins between and 6 GeV/c in 
the rapidity range \y\< 0.5. The combinatorial background was subtracted using like-sign 7T ± K ± pairs. 
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Fig. 3: (Colour online) (Upper panel) The invariant mass distribution in |y|<0.5 for the bin 0.4< 

/?T<0.5 GeV/c (left) and 0.9</?x<10 GeV/c (right), in pp collisions at 7 TeV. The background shape estimated 
using unlike-sign pairs from different events (event mixing) and like-sign pairs from the same event are shown as 
open red squares and full green squares, respectively. (Lower panel) The 7T ± K =F invariant mass distribution after 
like-sign background subtraction for 0.4</?x<0.5 GeV/c (left) and 0.9</?x<10 GeV/c (right). The solid curve is 
the result of the fit by Eq. Q] the dashed line describes the residual background. 
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Fig. 4: (Colour online) (Upper panel) The K + K~ invariant mass distribution in |y|<0.5, for the bin 
0.5</?t<0-6 GeV/c (left) and \.\<pi<\.2 GeV/c (right) in pp collisions at 7 TeV. The solid curve is the fit 
result (Eq.|2), while the dashed line describes the background. The background shape estimated using unlike-sign 
pairs from different events (event mixing) or like-sign pairs from the same event are shown as open red squares 
and full green squares, respectively. (Lower panel) The K + K~ invariant mass distribution after mixed-event back- 
ground subtraction for 0.5</?t<0.6 GeV/c (left) and l.l<px<l-2 GeV/c (right). The solid curve is the fit result 
(Eq. 13, while the dashed line describes the residual background. 
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Fig. 5: The product of acceptance and efficiency of K* and (j) detection as a function of pj in |y| <0.5. Statistical 
uncertainties are reported. Contributions to the point-to-point systematic uncertainties are listed in Tab. Q] 

The mass distribution M(n ± ,K T ) (see Fig. [3] for two pj bins) was fitted with a relativistic Breit-Wigner 
function multiplied by a Boltzmann factor Q and added to a polynomial residual background. The 
width was found to be compatible, within uncertainties, with the natural value. At low pj the fitted mass 
values were found to be slightly lower (by about «5 MeV/c 2 ) than the natural value, which is attributed 
to imperfections in corrections for the energy loss in the detector material. To extract the yield the 
distribution of M(n ± ,K T ) was then fitted with a (non-relativistic) Breit-Wigner function with the width 
fixed to its natural value (r= 48.7 ±0.8 fl25l ) and a background function: 



cW__l AT 

AM ~ In (M-M ) 2 + r 2 /4 + ( } ( ) 

where A is the area under the peak corresponding to the number of K* mesons, T is the full width at half 
maximum of the peak, and Mo is the resonance mass. The residual background B(M), after like-sign 
subtraction, was parametrized by a polynomial (dashed line in Fig. [3]). 

For the (j) meson, the raw yield was extracted from the K + K~ invariant mass distributions in 26 pj bins 
between 0.4 and 6 GeV/c. The combinatorial background was subtracted using a polynomial fit (first 
or second order), like-sign pairs, or unlike-sign pairs from mixed events (Fig. |4] for two pj bins). Since 
the invariant mass resolution of the <p peak is of the same order of magnitude as the natural <p width 
(~1 MeV/c 2 vs. 4.26 MeV/c 2 ), the fit is performed, after background subtraction, using a Voigtian 
function (convolution of Breit-Wigner function and Gaussian) superimposed on a polynomial to describe 
the residual background: 



— =A 't 9/ = AM' + B(M) (2) 

dM J {M-M') 2 + T 2 /A ^2Ha K ' 

where a represents the mass resolution and the other parameters have the same meaning as in Eq. £[]). 
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Fig. 6: Transverse momentum spectra for K* and 0(1020) in pp collisions at ^/s=l TeV. The statistical and sys- 
tematic uncertainties are added in quadrature and the uncertainty due to normalization |20] is shown separately. 
The statistical uncertainty is smaller than the symbol size. Each spectrum is fitted with a Levy-Tsallis function 
(dashed line). 



The background B(M) is represented in the lower panels of Fig. [4] by a dashed line. The width T is 
fixed to its nominal value [25] while a is a free parameter. The fitted mass values were found to be 
compatible, within uncertainties, with the known mass (25\ . with the exception of the low pi range 0.4- 
0.7 GeV/c where a fitted value lower than the natural one (by < 0.1%) was observed. The raw yields 
extracted using the three different methods to estimate the combinatorial background (analytic function, 
like-sign and mixed-event method) were found to be compatible within a few percent; therefore the mean 
value of all three methods was taken in each pj bin. 



4.4 Efficiency corrections 

In order to extract the meson yields, the raw counts (N RAW ) were corrected for the decay branching 
ratio [25 ] and for losses due to pion/kaon in-flight decays, geometrical acceptance, and detector efficiency 
(Af cor = N RAW /(A x e)BR, where BR indicates the decay branching ratio). The product of acceptance 
and efficiency (A x e) was determined for K* and <p from Monte Carlo simulations with the PYTHIA 6.4 
event generator (tune Perugia [14 ]) and a GEANT3-based simulation of the ALICE detector response. 
About 60 M Monte Carlo events, with the same vertex distribution as the measured events, were analyzed 
in exactly the same way as the data. The dependence on the event generator was estimated to be below 1 % 
by comparing PYTHIA and PHOJET simulated events. The Axe was determined from the Monte Carlo 
simulations as the ratio of the number of reconstructed resonances to the number of those generated, 
differentially as a function of rapidity and transverse momentum. The transverse momentum dependence 
is shown in Fig. [5] for K* and mesons. The decrease in A x e at low pj is due to the minimum 
pj requirement for reconstructed tracks, while the different behaviour for and K* is due to the different 
g-value of their decay (31.1 MeV for and 262.7 MeV for K*). 

Finally, corrections for the trigger efficiency (trigger) and the required primary vertex range (£ V ert) were 
applied in order to obtain the absolute resonance yields per inelastic collision: 
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Fig. 7: Energy dependence of (pj) for K* (triangles) and (j) (squares) in pp collisions. The points at lower energies 
are from STAR and PHENIX (^=200 GeV) EEES), ALICE (V?=0.9 TeV) [32) and E735 G/s=L8 TeV) l33l . 
The STAR data have been slightly displaced to separate the K* and the (p. The data point at 1.8 TeV represents the 
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d 2 N N™(pj) :: 1 x 
dydpx AyAp T e vert N M b 

here Af cor and Nmb are the number of reconstructed K* or <p and the total number of minimum bias trig- 
gers, respectively. The value of the trigger selection efficiency for inelastic collisions cogger is reported 
in Section 3. The loss of resonances due to the trigger selection, estimated by Monte Carlo, is negligible, 
less than 0.2%. The £ ver t correction factor accounts for resonance losses (s=T%) due to the requirement 
to have a vertex in the range of ±10 cm. 

4.5 Estimation of the systematic uncertainties 

The minimum and maximum values of the major contributions to the point-to-point systematic uncer- 
tainties are listed in Tab. [j] The uncertainty due to the raw yield extraction method was found to be 
±2-28% (2-10%) for K* (0). It was estimated by changing the mass range considered for the fit and 
the order of the polynomial for the residual background function (from first through third (second) order 
for K* (0)). Finally, variations in the yield due to the method used to estimate the combinatorial back- 
ground (like-sign and event-mixing method and also analytic function for 0) were incorporated into the 
systematic uncertainties. For the K* a relativistic Breit-Wigner function was used to fit the mass peak in 
addition to the non-relativistic version. In the case of the K* a rather large systematic uncertainty was 
estimated for the higher pj bins, due to the presence of a correlated background. 

The uncertainty introduced by the tracking and PID efficiency was estimated to be ±8% (8%) and ±1-6% 
(1.5%) respectively in the case of K* (0) by varying the kinematical and PID cuts on the daughter tracks. 
An additional ±1-4% uncertainty was added for the K* due to differences observed in the TOF matching 
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Table 1: Summary of the systematic point-to-point uncertainties in the K* and <j) yield 



Source of uncertainty K* 

Signal extraction ± 2-28 % ± 2-10% 

Tracking efficiency ± 8 % ± 8 % 

PID efficiency ±1-6% ±1.5% 

TOF matching efficiency ± 1-4 % 
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Fig. 8: (Colour online) Comparison of the K* pj spectra in inelastic pp collisions with PHOJET and PYTHIA 
tunes D6T (109), ATLAS-CSC (306), Perugia (320), and Perugia 201 1 (350). 



efficiency between data and Monte Carlo. The uncertainty on the yield contained in the extrapolated 
part of the <p spectrum was estimated to be ±20% using different fit functions. The normalization to the 
number of inelastic collisions leads to a +6% and -3% uncertainty in the yield of the measured particles. 
The resulting overall systematic uncertainty is tf% (t\ 2 { 5 %) for the K* (<j>) yield dN/dy and ±2% (3%) 
for the average transverse momentum {pj}- 
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Fig. 9: (Colour online) Comparison of the ( 1 020) pj spectra in inelastic pp collisions with PHOJET and P YTHIA 
tunes D6T (109), ATLAS-CSC (306), Perugia (320), and Perugia 201 1 (350). 
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Table 2: Parameters extracted from the Levy-Tsallis (O fits to the K* and (j) transverse momentum spectra in 
7 TeV pp collisions, including point-to-point systematic uncertainties. The first uncertainty is statistical and the 
second is systematic. 



Particles 


% 2 /ndf 


T (MeV) 


n 


K* 


2.0/19 


254±2±18 


6.2 ±0.07 ±0.8 





2.8/23 


272±4±11 


6.7 ±0.20 ±0.4 



Table 3: K* and <j> yield and (pj) estimated in the range 0-6 GeV/c in inelastic pp collisions at y/s=7 TeV. The 
systematic uncertainties of dN/dy and (pj) include contributions from the choice of spectrum fit function for 
extrapolation, the absolute normalization, and the point-to-point uncertainties listed in TabQ] 



Particles 


measured pj (GeV/c) 


dN/dy 


{p r ) (GeV/c) 


K* 



[0.0-6.0] 
[0.4 - 6.0] 


0.097 ±0.0004+g-g$ 
0.032 ±0.0004l°;j™ 


1.01 ±0.003 ±0.02 
1.07 ± 0.005 ±0.03 



5 Results and discussion 

5.1 px spectra and integrated yield 

Figure [6] presents the corrected pi spectra for the two resonances. The statistical and point-to-point 
systematic uncertainties added in quadrature are shown. The spectra are fitted with a Levy-Tsallis func- 
tion |[26ll27ll 

d 2 N (n-\)(n-2) dN ( m T -mV" 

— — = \ n ^tt x — x p T x 1 + — (4) 

dydpj nT[nT + m(n — 2 J J dy \ nT J 

where mj = \Jm 2 + p\. This function describes both the exponential shape of the spectrum at low py and 
the power-law distribution at large pj, quantified by the inverse slope parameter T and the exponent 
parameter n, respectively. The extracted parameter values are listed in Table [2] and the fits are shown 
in Fig. [6] The # 2 /ndf values are smaller than unity because the point-to-point systematic uncertainties, 
which are included in the fit, could be correlated. 

The extracted n values are similar to those quoted by the STAR experiment at RHIC for the measured 
in pp collisions at 200 GeV («=8.3±1.2) 0. In contrast, the slope parameters are significantly higher 
than the values obtained at RHIC, T= 202 ± 14 ± 11 MeV for 0, and T=223 ±8±9 MeV for K* (the 
latter was obtained by an exponential fit and can therefore not be directly compared). 

The total yields dN/dy and the mean transverse momentum (pj), including statistical and systematic 
uncertainties, are listed in Table [3] The values of dN/dy were obtained by integrating the spectra in the 
measured range and extrapolating to zero pj with the fitted Levy-Tsallis function. The contribution of 
the \ow-pj extrapolation is negligible for the K* and about 15 ± 3% for the (p. The mean transverse 
momentum was estimated in the range 0< pj < 6 GeV/c using the Levy-Tsallis function. However, 
similar values are obtained when calculating the mean from the measured data points, using the fit only 
to extrapolate into the unmeasured pj regions. In addition to the point to point systematic uncertainties 
previously described, an exponential fit was also used to estimate the systematic uncertainty in (pj) due 
to a different choice of fit function. Compared to pp collisions at 200 GeV |3j[5]|28]], the mean pj rises 
by about 30% (Fig. |7]) and the yield per inelastic collision increases by about a factor of two, which is 
similar to the overall increase of charged particle multiplicity |29l[30l . 
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Fig. 10: (Colour online) Energy dependence of the K*/K (upper panel) and <p/K* (lower panel) ratio in e + e 
(diamonds) J2]|4]]|42]|43|60]], and pp (triangles) fl] [3] |5] |47] |48] @9) collisions. Red squares represent the data 
from the ALICE experiment for 7 TeV pp collisions, K yields are from ifTTI . Open circles represent the same 
ratios in central nucleus-nucleus collisions from l3l 151 |49l l50l ISTI . Some points have been displaced horizontally 
for better visibility. Ratios are calculated from yields at mid-rapidity or in full phase space. 



The (j) yield, measured via the leptonic decay channel in the ALICE muon spectrometer in 2.5< y <4, 
l<px<5 GeV/c PTll . has a similar momentum distribution, but is lower by about 30% at forward ra- 
pidity. The <p yield is expected to vary by 20%-50% between forward (2.5 < y <4) and mid-central 
(-0.5< y <0.5) rapidities, based on analysis of different PYTHIA tunes described in paragraph 15.21 In 
particular, the lower value is predicted from the D6T PYTHIA tune 051 . which reproduces rather well 
the (j> spectrum at forward rapidity 11311 and the low pj part of the (j) spectrum at mid-rapidity (see Fig. [9] 
described in !5.2l) . 
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Fig. 11: (Colour online) Energy dependence of the <j) j% (upper panel) and (j) /K ratio (bottom panel) in nuclear 
(open circles) 0|2|5O]|52|53]|52|55], e+e" (diamonds) [EJ SH S3 Hgl [60]| , andpp (triangles) DJE|5]|28]|45]|47] 
E2|53] |54] ED collisions. Other n~ and K" yields are from fl56] E2J |58] [59J . Red squares represent the ALICE 
data at 0.9 and 7 TeV. The pion and kaon yields at 7 TeV are from ATI . The 0, 7ZT , and K~ yields at 0.9 TeV are 
from [24, 32 1. Some points have been displaced horizontally for better visibility. Ratios are calculated from yields 
at mid-rapidity or in full phase space, except the data at y/s=4. 87 GeV l52l . 



16 



The ALICE Collaboration 



5.2 Comparisons to models 

Multiparticle production, which is predominantly a soft, non-pertubative process, is usually modelled 
by QCD inspired Monte Carlo event generators like PHOJET @] and PYTHIA IM- In both models, 
hadronization is simulated using the Lund string fragmentation model [34]. Different PYTHIA tunes 
were obtained by adjusting the model parameters to reproduce existing data. The D6T tune ll35l . which 
uses the CTEQ6L parton distribution function (with a corresponding larger production of strange parti- 
cles), was obtained by fitting CDF Run 2 data. The ATLAS-CSC [36] tune was adjusted to minimum 
bias data from the UA5, E735, and CDF experiments for energies ranging from 0.2 to 1.8 TeV. The latest 
PYTHIA tune, Perugia 2011 lfT4l . takes into account first results from the LHC, in particular minimum- 
bias and underlying event data at 0.9 and 7 TeV. Strange baryon production was increased in this tune 
leading to a larger A/K ratio with respect to the Perugia tune. 

The transverse momentum spectra of K* and (j) are compared to PHOJET and various PYTHIA tunes in 
Figs. [8] and [9] For PYTHIA, tunes D6T (109), ATLAS-CSC (306), Perugia (320) and Perugia 2011 
(350) were used. The best agreement is found for the PYTHIA Perugia 201 1 tune, which reproduces both 
the K* spectrum and the high pi part (pr> 3 GeV/c) of the <p spectrum rather well. PHOJET and ATLAS- 
CSC very significantly overestimate the low momentum part (pj< 1 GeV/c) of the transverse momentum 
distribution but reproduce the high momentum distribution of both mesons well. The PYTHIA D6T tune 
gives the best description at low pj, but deviates from the data at pi>2 GeV/c. Finally, the PYTHIA 
Perugia tune underestimates the meson yield for pj larger than 0.5 GeV/c. 

Similar comparisons for the mid- and forward-rapidity spectrum in pp collisions at a/s=0.9 TeV 11321 
and 7 TeV OTI . respectively, show that the spectrum is rather well reproduced by the ATLAS-CSC 
and D6T tunes, while the Perugia and 201 1 tunes underestimate the data. Moreover the PYTHIA tunes 
generally underestimate strange meson and hyperon production in 7 TeV pp collisions Ifl2"ll37l . while 
the Perugia 201 1 tune gives a good description of kaon production in pp collisions at 7 TeV [11]. 

5.3 Particle ratios 

The measurement of particle production and particle ratios in pp collisions is important as a baseline 
for comparison with heavy ion reactions. In heavy ion collisions, the yields for stable and long-lived 
hadrons reflect the thermodynamic conditions (temperature, chemical potentials) at freeze-out, whereas 
the yield for short-lived resonances can be modified by final state interactions inside the hot and dense 
reaction zone 1381 1391 . Particularly interesting is the comparison of and K* production, considering 
the different lifetimes (about a factor 10) of the two resonances. 

Using different particle ratios (like K/n or 0/K*) measured in elementary collisions, values ranging 
from 0. 1 to 0.4 fl] |2j |41] |42j |43j were previously obtained for the strange quark suppression factor 
A s = 2ss/(uu + dd), which represents the probability to produce strange quark pairs relative to light 
quarks BOl . In pp reactions, particle abundances have been successfully described by statistical-ther- 
mal models. Now, using measured identified particle yields, an energy-independent value of 0.2 for A s 
has been extracted in e + e~, pp, and pp collisions at y / ?<l TeV Il40ll44l . 

Using the <p and K* yields presented in this paper and stable particle results measured by ALICE at the 
same energy ifTTTl . we find the following values for particle ratios in pp collisions at 7 TeV: K*/K~= 
0.35± 0.001 (stat.) ±0.04 (syst.), <j)/K*= 0.33 ±0.004 (stat.) ±0.05 (syst), 0/Kr= 0.11 ±0.001(stat) 
±0.02 (syst.), 0/7T~= 0.014 ±0.0002 (stat.) ±0.002 (syst.). Due to the fact that the same data were 
analyzed to extract both resonance and non-resonance (jt,K) yields, the uncertainties due to the absolute 
normalization cancel and are therefore not included in the systematic uncertainties of the ratios. These 
ratios are shown in Figures [10] and [TT] together with the results obtained at lower incident energies in pp, 
e + e~, and A-A collisions. 

The K*/K~, 0/K~, and (j>/K* ratios are essentially independent of energy and also independent of the 
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Fig. 12: (Colour online) (£2 + + Q. )/0 ratio as a function of transverse momentum for pp collisions at y/s=7 TeV. 
D. data are from [ 12 1. The dashed line represents the prediction of HIJING/BB v2.0 model with a SCF for pp col- 
lisions at -^=5.5 TeV with a string tension of 2 GeV/fm [13|. The same calculation at 7 TeV yields a ~ 10% 
higher ratio [62 1 . The full line represents the prediction of the PYTHIA Perugia 201 1 tune [14| for pp collisions at 
V^=7 TeV. 



18 



The ALICE Collaboration 



collision system, with the exception of K7K and 0/K* at RHIC (5] |49j [50l 152, where these ratios in 
nuclear collisions are respectively lower and higher than in pp. On the contrary, the /% ratio increases 
with energy both in heavy ion and in pp collisions up to at least 200 GeV. However, in heavy ion collisions 
the value obtained by the PHENIX experiment B, about 40% lower than the STAR result [5] at the 
same collision energy, seems indicate a saturation of this ratio at the RHIC energies. In pp collisions we 
observe a saturation of the <p /% ratio, with no significant change over the LHC energy range between 1 
and 7 TeV. 

In microscopic models where soft particle production is governed by string fragmentation, strange hadron 
yields are predicted to depend on the string tension [8]. Multi-strange baryons, and in particular the ratio 
are expected to be very sensitive to this effect lfT3l . The yield is compared to the D. data 
measured by ALICE at the same incident energy [12] in Fig. [12] as a function of transverse momentum. 
The full line represents the PYTHIA model (Perugia 2011 tune), which is a factor 1.5-5 below the data. 
While this tune describes the <p spectrum reasonably well above 2-3 GeV/c, it underpredicts multistrange 
baryon yields by a large factor [12]. The dashed line, which is very close to the data, represents the 
prediction of a model with increased string tension, the HIJING/BB v2.0 model with a Strong Colour 
Field (SCF), for pp collisions at 5.5 TeV lPT3l . This is a model that combines multiple minijet production 
via perturbative QCD with soft longitudinal string excitation and hadronization. In this case the SCF 
effects are modeled by varying the effective string tensions that controls the qq and qqqq pair creation 
rates and the strangeness suppression factor. The value of string tension used in this calculation is fc=2 
GeV/fm, equal to the value used to fit the high baryon/meson ratio at -^=1.8 TeV reported by the CDF 
collaboration [61]. The same calculation at 7 TeV yields a ~ 10% higher ratio (62J. Higher values of the 
string tension (~ 3 GeV/fm) also successfully reproduce also the (Q. + £2) /<j) ratio in Au-Au collisions at 
V^=200 GeV (131, but overestimate the (A+ A)/K| at 7 TeV 

6 Conclusion 

Yields and spectra of K*(892)° and 0(1020) mesons were measured for inelastic pp collisions at y/s=l TeV 
by the ALICE collaboration at the LHC. The transverse momentum spectra are well described by the 
Levy-Tsallis function. The yields for both mesons increase by about a factor of two from 200 GeV 
centre-of-mass energy, and the average pj by about 30%. 

The K*/K and 0/K* ratios (and consequently the 0/K ratio) are found to be independent of energy up to 

7 TeV. Also the (p/n ratio, which increases in both pp and A- A collisions up to at least RHIC energies, 
saturates and becomes independent of energy above 200 GeV. 

The data have been compared to a number of PYTHIA tunes and the PHOJET event generator. None of 
them gives a fully satisfactory description of the data. The latest PYTHIA version (Perugia 2011) comes 
closest, while still underpredicting the <p meson pj spectrum below 3 GeV/c by up to a factor of two. 

The (Q.~ + Q + )/(j) ratio is not reproduced by PYTHIA Perugia 2011, but is in good agreement with the 
HIJING/BB v2.0 model with SCF, which enhances multi-strange baryon production by increasing the 
string tension parameter. 
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